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VQ , The structure of a molecular liquid, in both the nematic liquid crystalline and isotropic phases, 

f^ ' around a cylindrical macroparticle, is studied using density functional theory. In the nematic phase 

C^ ' the structure of the fluid is highly anisotropic with respect to the director, in agreement with results 

^S| . from simulation and phenomenological theories. On going into the isotropic phase the structure 

, , ■ becomes rotationally invariant around the macroparticle with an oriented layer at the surface. 
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PACS numbers: 61.20.Gy,61.30.Cz,61.30.Jf 



I. INTRODUCTION 

c/j ■ Solid particles, both spherical and non-spherical, dispersed in a liquid crystal (LC) host comprise an interesting 

class of novel materials [llQl- In orientationally ordered phases of the host, the introduction of solid particles deforms 
Ci I the director field, leading to long-range interactions between the particles and effects such as chaining or the formation 
of soft solids. Colloidal dispersions in liquid crystals have a wide range of applications 3] and have recently attracted 
a great deal of interest ,Jj . 

Experimental techniques such as atomic force microscopy or confocal microscopy may be used to study LC-colloid 
dispersions; simulations and theory have also been applied. Simulations have been used to study the ordering of LC 
fj molecules around one ^, ^J or two ma crop articles |3 , and LC dispersions in confined geometries ^] . Phenomenological 
theories such as Landau-de Gennes 1^113 Ell IQ O'^ Frank elastic [13,13 theory have also been used. These approaches 
however have their limitations: simulation is computationally expensive, while the aforementioned phenomenological 
theories require, often poorly known, parameters, and are incapable of accounting for spatial variation in the density 
tH- ■ (and, in the case of elastic theory, variation in the order parameter). 
T^ ' One popular theoretical method that may be applied to this type of problem is density functional theory (DPT) 

C^l . UM- Unlike Landau-de Gennes or clastic theory, this is capable of accounting for spatial variation in the density 
\l ■ and order parameter and, in the form used here, requires knowledge only of the interaction potential between the 
C^ . molecules in the fluid. In this case DPT at the level of the Onsager second virial approximation [l^ is applied to 
the case of a single, infinitely long, cylindrical macroparticle in a LC host. As the system is homogeneous along the 
length of the cylinder, two dimensions are sufficient to represent the spatial dependence of the density. 



II. THEORY 

^ ■ 

C . Por a system of uniaxial molecules the grand potential can be written as [ij 

O 



(3n [p{r, u)] = (3Fi^ [p{r, u)] + (3F,^ [p(r, u)] + f3 drdu {V,^t{r, u) ~ p)p{r, u) , (1) 



where p{r,u) is the position- and orientation-dependent single particle density, Vcxtir,u) is the external potential, 
j^ ' /i is the chemical potential, and (3 = 1/kBT. F\^[p{r,u)] and Fcx[p{r,u)] are the ideal and excess free energies, 
respectively, r is the position vector and u is the orientation vector. The ideal free energy is given by 

(3Fid[p{r,u)]= drdu p{r,u) {log p{r,u)-l} . (2) 
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The exact form of the excess free energy is generally unknown. Here we employ the Onsager approximation [13 

/5Fex[p(r,M)] = -- / dridr2duidM2 /(ri2,Mi,M2)p(ri,Mi)p(r2,M2) , (3) 

where /(ri2,Mi,M2) = exp{— /3y(ri2, Mi,M2)} — 1, the Mayer function, ri2 — ri ~ r2 and y(ri2, 1^1,1*2) is the 
molecular pair interaction potential. The liquid crystal here is modelled as a fluid of prolate hard ellipsoids of 
elongation e — a/b = 5; a is the length of the symmetry, or major, axis and we will use the minor axis 5 = 1 as a 
unit of length. When two molecules overlap, F = 00 and / = —1; for a non-overlapping pair, V = and f — 0. The 
approximation of eqn ^ corresponds to truncating the virial expansion after the pair term. While Onsager theory is 
exact only in the limit of infinite elongation, it has been used to study the anchoring of ellipsoids of this elongation 
near solid substrates p^ll9| and has been found to give results in qualitative agreement with simulation. 

The external potential, representing a single cylindrical macroparticle of radius R oriented along the y axis, is given 
by 



Ve^t{r,u) = Ve^tis,u) = I Wo 



\Vq [tanh(6/-u;) - tanh(-6/u')] s - R < -b 

tanhf ) — tanh(— 5/u;) \s 

\ w / 



-R\<b (4) 



w 
s-R>b 

where s — (x, z), s — \s\, Vq — SO/cbT and w — 6/5. This represents a sharply varying repulsive potential acting 
on the ellipsoid centres of mass. It excludes the molecules from the cylinder and gives rise to homeotropic (normal) 
anchoring at the surface. 

As before j20. 21 . 22. 23] the angularly dependent functions are expanded in a set of spherical harmonics, and the 
assumption of translational invariance along y allows us to write the coefficients as functions of s: 

log p{r,u)^^Pijn{s)Yem{u) (5a) 

i,m 

P{r,u) = ^p^™(s)Y^;„(m) (5b) 

l/oxt(r, m) = ^ Vtm{s)Yi„i{u) . (5c) 

Note the complex conjugate in the density expansion. The Mayer function is expanded as |22| 

/(ri2,Ml,M2) = ^ flil2l!.{'ri2)'^lit2t{'ri2-,Ui,U2) , (6) 

where ri2 = |^i2|, 1^12 = ^i2/''i2, and ^i^e,2t is a rotational invariant J2J|. Inserting these expressions into the grand 
potential and integrating over angles and the y direction gives 



l3Vt[p{r,u)] 



/ ds ^ Pf,,n{s) [ptm{s) - V47r(l + Pp)5l0 + (3Vlra{s) 
/ dSidS2 ^ Aimi£2m2(si2)P£imi(si)/3f2ni2(s2) ■ 



(7) 






Here L is the box length in the y direction (we assume periodicity). The quantities Ci^rnxi^ra^i^vi) come from 
integrating the Mayer function and are the spherical harmonic coefficients of the excluded length (in the y direction) 
of two molecules with a separation vector S\2 = Si — S2 in the xz-plane, treated as a function of the molecular 
orientations. As the last term in eqn Q is a convolution, it is most conveniently evaluated in reciprocal space. If 
Pim{k) is the two-dimensional Fourier transform of pimis) then this term may be written 



^ ^ Ci^mie2ni2{l^)Peinii{k)pi2ni2{k) (8) 

k iimi 

^21^12 

where Ci^mii2m2{k) is the Fourier transform oi Ct^rnii2m2{^i2)- 



In order to find the equilibrium density the functions are tabulated on a regular grid in the xz plane; the grid spacing 
is 6x ~ Sz ^ 0.26, the molecular length corresponding to 25 grid points. The grand potential is then minimised with 
respect to the pimis) coefficients at each grid point using the conjugate gradient method J23. When required, the 
coefficients phn{s) are calculated through eqns jSJl, with angular integrations performed using Lebedev quadrature 
|26l l27l|. The systems in studied in this work were all square boxes {L^ — Lz) with L^ ranging from 406 (200 grid 
points) to 1006 (500 grid points). Table Q] summarises the parameters of the systems discussed in detail in the next 
section. 



R/h Lx/b Ux 8kx/h ^ 
5 40 200 7r/20 
7.5 60 300 7r/30 
10 60 300 7r/30 
15 80 400 7r/40 
20 100 500 7r/50 

TABLE I: Details of systems studied. R is the macroparticle radius, L^ is the box length in the x direction (with L^ = L^), 
Ux is the number of grid points in the x direction and 5kx is the grid spacing in k space. 

Once the equilibrium density coefficients pim{s) have been determined, the number density p[s) around the 
macroparticle may be found from 

p[s) = / dit p{s,u) = %/47rpoo(s). (9) 

The orientational ordering is described by the order tensor Qa(j{s) that is found from p{s,u) by 

Qapis) = ^ du p{s,u)ua{s)up{s) - ^Sa0, a,(3 = x,y,z. (10) 



The spatially varying order parameter S{s) is given by the largest eigenvalue of Qaf3{s) and the director n(s) by the 
eigenvector associated with S{s). 

III. RESULTS 
A. Structure in the nematic phase 

First we examine the fluid structure around the cylindrical particle in a nematic fluid, at chemical potential p = 8.0. 
The density distributions around cylinders of radius 5 < R/b < 20 are shown in Fig. ^ For all R the density is largest 
at the surface but then decays away with almost periodic variation, similar to that seen for a nematic-planar wall 
interface. Further from the cylinder this distortion in the density becomes highly anisotropic. Parallel to the director 
the modulations in the density are stronger than in the perpendicular direction. This weakening perpendicular to the 
director is due to the partial melting of the nematic in the defect regions [23 and is most noticeable for the smallest 
radius R/b = 5 (Fig. ^). In this case, there is almost no density modulation perpendicular to the director causing 
the density map to show chevron-like structures, with the chevron tips pointing along the director away from the axis 
of the cylinder. For larger radii, the density variation perpendicular to the director is stronger, but shorter-ranged, 
than in the direction parallel to the director. In comparison to simulation of a system with e = 3 |28J, the density 
variation seems to be shorter ranged. While this may be due to differences in the models, studies of the nematic-wall 
interface using Onsager and related theories [i&L |2^ also gave density profiles that have generally weaker structure 
than comparable simulations. 

Maps of the orientational order parameter around the cylindrical particles are shown in Fig. |21 As for the density 
maps, the variation in the order parameter is highly anisotropic. For R/b = 5, along the director, there are lobes 
of high order along the top and bottom of the cylinder. In the direction perpendicular to the director there are two 
regions of drastically reduced order, corresponding to defects in the liquid. As R increases the lobes of increased 
order tend to wrap around the cylinder, and the defects move away from the surface, in agreement with simulation 
[23 and phenomenological theory ^^ l29l 1301 . While the positions of the defects are in qualitative agreement with 
previous results, in the present case the defects are significantly smaller than those seen previously, suggesting that for 
molecules of this elongation more sophisticated methods such as weighted density approximations 31] or fundamental 
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FIG. 1; Density maps around a cylindrical macroparticle at /i = 8.0 (neniatic phase). Cylinder radius Rjh: (a) 5.0, (b) 7.5, (c) 
10.0, (d) 15.0, and (e) 20.0. Dark colours show areas of low p(s), light colours high pi^s). 
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FIG. 2: Order parameter maps around a cylindrical macroparticle at ^ = 8.0 (nematic phase). Cylinder radius R/b: (a) 5.0, 
(b) 7.5, (c) 10.0, (d) 15.0, and (e) 20.0. Dark colours show areas of low S{s), light colours high S{s). Solid lines show the 
orientation of the local director n{s). 



measure theory [33 are necessary to examine the structure within the topological defect. Detailed comparisons with 
simulation will appear in a later publication. 

Also shown in Fig.|51is the director orientation n(s) around the macroparticle. At large distances from the cylinder 
the n{s) lies along the z axis. Close to the particle the director becomes highly distorted. Along the z axis and at 
the defects n{s) is normal to the particle surface. At other points on the particle surface, n(s) points away from 
the surface normal appearing to graze the surface of the macroparticle. This behaviour is different from that seen in 
simulation z&| and from elastic theory in the case of strong anchoring 29] , suggesting that the anchoring around the 
macroparticle is weaker than in previous studies. Apart from in the vicinity of the defects n(s) is a smoothly varying 
function of s. 

In order to gain more insight, the density and order parameter profiles parallel (py (s) and Sh (s)) and perpendicular 
{p±{s) and S±{s)) to the director (through the defects) are shown in Fig.|31 The p\\{s) curves show similar structure 
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FIG. 3; 

(a) Density profiles, p{s)b , and (b) order parameter profiles S{s), around a cylindrical macroparticle at ^ = 8.0 (nematic 

phase). Full lines: in the direction parallel to the director. Dashed lines: in the direction perpendicular to the director. 

Results are shown for cylinder radius R/b = 5.0 and, displaced successively upwards by 0.1 units for clarity, 

R/b^ 7.5,10.0,15.0,20.0. 

as R increases, with a peak at the cylinder surface s = R, a second peak about 56 (one molecular length) from the 
surface, and some decaying oscillations into the bulk. This is similar to the behaviour seen for a nematic fluid near 
a wall with homeotropic alignment 20, 23]. The size of the peak at contact increases with R at small R. p±{s) also 
has a peak at contact. For smaller cylinders p±{R) <C P\\{R)- As R increases the difference between the two heights 
decreases. Further from contact pj_{s) has a secondary peak. For the R/b — 5 cylinder this peak is approximately 
2.56 from the surface with p±{s) tending towards the bulk value further out. For larger R, the secondary peak moves 
out. Beyond this peak p±{s) behaves almost identically for all R. 

As with /3||(s), the S'||(s) profile is similar to the order parameter profile for a nematic-homeotropic wall interface, 
with a maximum at contact and decaying oscillations. S±(r) shows a minimum near contact. This corresponds to the 
defect seen in the order parameter maps (Fig. [21 . In agreement with the order parameter maps, the distance between 
the minimum and the surface increases with R. Shown in Fig. ^ is the position of the minimum Smin as a function of 



radius R. For spherical macroparticles with both Saturn ring and satelHte defects, the distance between the surface 

T 




R/b 

FIG. 4: Defect position (smin — R) near cylindrical macroparticle at /i = 8.0 (nematic phase) as a function of radius R. Dashed 
line shows the line of best fit. 

and the defect varies linearly with particle size ^(f\ . In the present case s„iin increases with R. 

B. Structure in the isotropic phase and around the nematic-isotropic transition 

Shown in Figs. ^^, \2^, El andEI are the density and order parameter maps for a cylindrical macroparticle with 
R/b — 15 in both the nematic and isotropic phases. The bulk isotropic-nematic transition occurs at /i^v/ = 7.3675. 
In the nematic phase, at this value of /i, the density and order parameter maps are similar to those seen for /z = 8, 
though the variation in the density and order parameter are less pronounced. Also the defects either side of the 
cylinder with /i = 7.3675 are both larger and significantly further from the surface (smin — R ~ 2.5b) than for the 
higher chemical potential. In the isotropic phase, both the density and order parameter variation become almost 
completely symmetrical, reflecting the loss of a preferred direction. In this phase, far from the cylinder, the order 
parameter 5 — > 0, while both p and S have maxima at the surface. This suggests that the surface is wet by the 
nematic and this give s rise to short range interactions between particles dispersed in a liquid crystal host even in the 
isotropic phase pAl33|. However, a detailed investigation of surface phase behaviour and interparticle interactions in 
this model await further study. 



IV. CONCLUSIONS 



In this paper the structure of a liquid of hard ellipsoids of elongation e = 5, around a cylindrical macroparticle, in 
both the nematic and isotropic phases, was studied using density functional theory within the Onsager approximation. 
The resulting density and order parameter maps were consistent with previous theoretical and simulation work. On 
going from the nematic to isotropic phase, the structure of the surrounding fluid becomes rotationally invariant about 
the cylinder, with what appears to be a nematic wetting layer at the particle surface. 

The present study is preliminary: investigation of the sensitivity of the results to the resolution of the real-space 
and reciprocal-space grids must still be carried out. Nonetheless, the results are very promising. Despite its simplicity. 




FIG. 5: Density maps for fluid around cylindrical macroparticle of radius R/b — 15 with (a) /i = 7.3675 (nematic phase, at NT 
transition), (b) /i = 7.3675 (isotropic phase, at NI transition), and (c) jj, = 5.0 (isotropic phase). Dark colours show low p{s), 
light colours high p{s). 

this method provides results that are generally in agreement with simulation and phenomenological theory. The only 
major deficiency is that the defects are smaller than would be expected from simulations of a similar system with 
e = 3. Using more sophisticated density functionals or ellipsoids with longer elongations (where Onsager theory is 
more accurate) would hopefully give a better description of the defect. Detailed comparisons with simulation will be 
described elsewhere. 
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